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Abstract 

A  scalar  plane  wave  incident  on  a  penetrable  sphere  is  considered 
in  the  short  vave  length  limit.  A  new  representation  of  the  scattering 
amplitude  is  introduced  which  is  particularly  appropriate  in  this  limit, 
and  which  req\ilres  only  the  evaluation  of  certain  integrals.   Some  of 
these  may  be  evaluated  asymptotically  by  the  method  of  steepest  descent 
and  lead  to  the  geometrical  optics  field  contribution.   Included  in  this 
Is  the  bow  field.   The  remainder  of  the  integrals  are  evaluated  by  the 
method  of  residues  and  lead  to  the  diffracted  field  contribution.   This 
"diffracted  ray"  field  is  known  from  recent  investigations  in  diffrac- 
tion theory.   An  essential  part  of  the  analysis  is  the  introduction  of 
the  parameter  p,  the  number  of  internal  refractions  that  a  ray  which 
hits  the  sphere  undergoes.   The  results  obtained  ape  all  in  agreement 
with  that  which  would  be  expected  on  the  basis  of  geometrical  diffrac- 
tion theory. 
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I.   Introduction 

We  consider  a  scalar  plane  wave  incident  on  a  penetrable  sphere  in  the 
limit  of  high  frequencies  for  which  the  wave  length  of  the  incident  wave  is 
small  compared  to  the  radliis  of  the  sphere.   Quajitimi  mechanically,  the  prob- 
lem represents  the  scattering  of  a  highly  energetic  particle  by  an  attractive 
square  well.   It  is  also  closely  related  to  the  corresponding  electromagnetic 
problem  which  occurs  in  the  scattering  of  light  by  water  droplets.  As  such, 
it  has  engaged  the  attention  of  very  many  scientists  whom  we  can  not  hope  to 
enumerate  here.  A  more  or  less  complete  historical  record  of  the  contributors 
to  the  raindrop  problem  may  be  found  in  van  de  Hulst's  book  ^ -^  .      The  Jiistifi- 
cation  of  the  present  work  lies  in  the  new  representation  of  the  scattering 
amplitude  which  is  particularly  appropriate  to  the  study  of  it  in  the  high 
frequency  or  short  wave  length  limit.  By  means  of  it,  it  is  possible  to  dis- 
tinguish the  contributions  of  the  geometric  optics  field  and  of  the  diffracted 
field. 

In  a  previous  report  '-  -■  the  forward  scattering  amplitude  alone  was  consi- 
dered. In  the  present  study  we  consider  the  scattering  amplitude  at  arbitrary 
angles.   The  basic  idea  utilized  is  the  same  as  that  of  Rubinow  and  Wu  ^  -^    and 
Wu  i- -■  in  treating  the  forward  scattering  from  impenetrable  cylinders  and  spheres, 
It  consists  in  converting  the  Infinite  series  representation  of  the  scattering 
amplitude  into  a  sum  of  integrals.  We  then  utilize  an  expansion  first  intro- 
duced by  Debye  in  I908  '--'-'  in  his  work  on  the  dielectric  cylinder.  In  fact, 
it  is  worth  remarking  that  the  spirit  of  this  work  is  no  different  from  Debye 's 
work.   By  means  of  Debye' s  expansion,  it  is  possible  to  express  the  scattering 
amplitude  F(9)  a£  a  sum  over  p  of  partial  amplitudes  F  (9).   The  geometric  mean- 
ing of  p  is  that  it  represents  the  number  of  refractions  that  a  ray  which  hits 


the  sphere  undergoes.  The  partial  amplitude  F  (9)  represents  the  contribution 
to  the  scattering  at  angle  9  of  the  particular  ray  characterized  by  p  refrac- 
tions. 

Each  partial  amplitude  is  represented  by  a  sum  of  integrals.   It  is  then 
possible  to  rearrange  the  integrals  so  that  they  split  into  two  groups.   In  one 
group,  the  integrals  have  stationary  points  along  their  integration  paths.   In 
the  second  group,  they  do  not.   In  making  this  split  up,  an  idea  introduced  by 
Franz  ^  J  is  useful.   The  first  group  of  integrals  may  be  evaluated  asympto- 
tically by  the  stationary  phase  method  or  method  of  steepest  descent.   The  result 
is  called  the  geometric  optics  field  because  it  is  just  what  would  be  expected 
from  purely  geometric-optical  considerations.   The  second  group  of  integrals, 
called  the  diffracted  field,  is  evaluated  by  the  method  of  residues.   They 
correspond  physically  to  waves  or  rays  which  encircle  the  sphere  and  decay  expo- 
nentially.  Such  phenomena  are  known  from  the  work  on  impenetrable  cylinders 
and  spheres  of  Franz  and  Deppermann  L'-i,  where  they  are  called  creeping  waves  . 
They  are  also  known  from  the  geometrical  theory  of  diffraction  of  Keller  <-   -J, 
in  which  they  are  more  appropriately  named  "diffracted  rays".   They  occur  in 
the  scattering  from  smooth  three  dimensional  objects  L-^J  as  well  as  from  convex 
cylinders  L  -I  . 

In  this  paper  the  language  of  geometrical  diffraction  theory  is  used  through- 
out. Not  only  does  this  permit  a  ready  physical  identification  of  the  various 
formulas,  but  likewise  it  provides  a  convenient  physical  check  on  the  work.   The 
rainbow  phenomenon  appears  here  as  a  result  of  a  stationary  point  of  the  third 
order.  For  this  reason,  there  is  a  greater  intensity  of  the  field  at  the  bow 
angle.   The  bow  angle  is  also  a  caustic  of  the  p-family  of  rays,  and  is  a  local 
minimum  or  maximum  of  the  scattering  angle. 
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Yet  to  be  examined  is  the  field  in  the  neighborhood  of  the  axial  caustics, 
i.e.  in  the  forward  and  backward  directions,  as  well  as  the  field  in  the  neigh- 
borhood of  the  boimdary  of  the  illuminated  region  associated  with  a  given  value 
of  p. 

II.   Formulation 

A  scalar  plane  wave  traveling  in  the  direction  9  =  0  is   incident  on  a 
penetrable  sphere  of  radius   a.      The  wave   function  \\i  satisfies   the  scalar  helm- 
holtz  equation 

A^+  n\\  =  0  r  ^   a    ^ 

''' 

A\)/+  kt=0  rSaj. 

The  quantity  N  is  called  the  refractive  index  of  the  sphere.   The  wave  function 

and  its  normal  derivative  must  be  continuoiis  at  r  =  a.  At  large  distances  from 

the  sphere,  the  wave  function  consists  of  the  incident  wave  plus  a  scattered 

field: 

^1      o        ikr 

ikr  cos  9   „/„N  e  ,  , 

t  -e         +  f(9)  -^~   ,  (2) 

The  function  of  interest  f(9)  is  called  the  scattering  amplitude.  As  a  con- 
venience, we  introduce  the  non-dimensional  scattering  amplitude 

F(9)  =  f(9)/a  (3) 

and  non-dimensional  parameter 

x  =  ka  .  (4) 

The  problem  as  defined  represents  the  quant-um-mechanical  non-relatlvistic 


scattering  of  a  particle  by  a  square-well  potential.      If  the  well  Is   attrac- 


N  =      \/l  +  2mV  a^  /  Il^x^ 


(5) 


with  m  the  mass  of  the  particle  and  h  equal  to  (ait)   times  Planck's  constant. 
Nevertheless,  we  find  it  more  descriptive  to  employ  the  language  of  geometric 
optics  and  diffraction  theory.   In  other  words,  we  think  of  this  problem  as 
representing  the  scattering  of  light  from  a  dielectric  sphere,  idealized  to 
the  case  in  which  the  Incident  wave  is  a  scalar.   For  simplicity,  N  will  be 
assumed  to  be  real  (i.e.,  non-absorbing  sphere)  and  greater  than  unity. 

The  scattering  amplitude  has  the  well  known  series  representation  ' ' 

F(9)  =  I  2  (n  +  l/2)P^(cos  9)  [l  -  e^^^^^^J  (6) 


vhere  ,    .  ,    ., 

^i26(n)  ^  _  ^r(-)JA(^)  -  y  /^)-^n(^) 
Nh^^^(x)j;(Nx)  -h^^^'(x)jjNx) 


(T) 


P  is  the  Legendre  polynomial  of  degree  n,  J  ,h  are  spherical  Bessel,  Hankel 
n  n  n 

functions,  respectively,  and  primes  denote  differentiation  with  respect  to 
argvment . 

We  seek  the  asymptotic  form  of  F(0)  for  large  values  of  x.   This  is  the 
neighborhood  of  the  geometric  optics  limit  for  which  the  wave  length  (2jt/k)  of 
the  incident  wave  is  small  compared  to  the  size  of  the  scattering  object.   To 
examine  this  behavior,  the  representation  (6)  is  not  at  all  useful,  because 
ver^'  many  terms  contribute  to  the  sum,  the  number  of  such  terms  being  of  the 
order  of  x.  We  utilize  Instead  a  representation  which  essentially  converts 


the  Slim  into  an  integral  over  the  index.  This  conversion  has  been  utilized 
successfully  by  Rubinow  and  Wu  '-^'-1  ^  and  Wu  L  J  ^  in  finding  the  correction  to 
the  geometrical-optics  scattering  from  impenetrable  cylinders  and  spheres. 

ri2i 

The  method  is  not  unrelated  to  the  method  of  the  Watson  transform^   ,  which 
has  also  been  used  extensively  in  recent  years  in  the  investigation  of  dif- 
fraction problems. 


Ill   The  Scattering  Amplitude  Representation 

We  introduce  the  Fourier  transform  of  the  summand  in  equation  (6).   The 
summand  is  assumed  to  be  defined  for  n  >  -  l/2.  By  means  of  the  Poisson  sum- 
mation formula  L  J^  equation  (6)  becomes 


m   =  i  Z^  /  ^v  V  P,.,/2(=-  «)[l  -  ,'^(^-^/^^y^'-(^-^/^K        (8) 


The  function  P    /2( cos  G)  is  to  be  understood  as  the  hyper geometric  function 
which  is  identical  with  the  Legendre  fimction  only  for  integral  values  of  v. 
It  is  singular  at  9  =  it.  Although  the  Poisson  formula  effects  the  conversion 
from  simmation  over  the  discrete  index  n  to  integration  over  the  continuous 
variable  v,    an  infinite  sum  remMns,  now  taken  over  the  new  index  m.  However, 
it  will  be  seen  that  the  major  contribution  arises  from  the  m  =  0  term,  just 
as  for  the  conducting  sphere  problem.   We  desire  to  evaluate  the  integral  for 
large  x  by  the  method  of  steepest  descent.   If  the  phase  function  6(v-l/2)  is 
examined  as  in  reference  [jj ,  it  will  be  found  that  the  Integrand  in  question 
has  very  many  stationary  points  lying  on  the  real  axis.   The  Integrand  can  be 
expanded  in  such  a  way  that  each  term  contains  only  one  stationary  point.   In 
view  of  the  relation 
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KM 


\l  2x. 


V+l/2 


(x) 


(9) 


for  Hankel  functions   of  either  the  first  or  second  kind,   ve  have 


(2) 
i2S(v-l/2)   _  \      M 


\^  +  \2  ^21     2     Rgf" 


LHi^^(Nx). 


(10) 


vhere 


^11  = 

U 
"  L 

^22  = 

V 
"  L 

^12  = 

1  + 

^1 

= 

^21  = 

1  + 

«22 

= 

h\ 


4i 
rtxHj,^^(Nx)H^^Vx)L 


.(2) 


(x) 


.(2)^ 


(Nx) 


'(x) 


J^ 


(Nx) 


h(^)'(x) 


(1) 


(Nx) 


-m 


(x) 


7^^ 


H^^^'(Nx) 


X2T 


(Nx) 


(11) 


This  expansion  vas  first  Introduced  by  Debye  L-^-l .  The  quantities  R  ^  Rp_  are 
intei-pretable  as  Fresnel  reflection  coefficients  for  the  interface  (r  =  a)  be- 
tween the  exterior  (subscript  l)    and  the  interior   (subscript   2)    of  the   sphere  L     •*    -^-' 
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Thus,  R   represents  the  reflection  coefficient  of  an  inward  spherical  wave 
incident  from  the  exterior  emd  reflected  back,  while  R^p  represents  the  re- 
flection coefficient  for  an  outward  spherical  wave  emanating  from  the  center 
of  the  sphere  and  reflected  back.   The  quantity  T^ _  is  the  transmission  coef- 
ficient for  the  wave  moving  from  the  exterior  to  the  interior,  and  Tp,  is  the 
transmission  coefficient  for  the  reverse  process. 

After  inserting  equation  (lO)  into  (8),  we  obtain 


F(Q)  =  S  F  (0)  .  (12) 

p=o 


2jtmv 

"v-1/2^ 


O  l_       y     \   '      _I 

00  00 

V^)  =  "  X  S^  f    ^'   ^p,m  Vl/2(^°^  ^)     P  =  1^2,5,...  .     (11. 

O 


(2),„.  n,(i)(,,p^ 


^p,.=  v(-)"  T,,  ^r  -..  irn-     iar 


(Nx) 


12nmv 
p  =  1,2,5,...  .     (15) 


It  is  seen  that  F(9)  decomposes  into  parts  labeled  Fq(9)  and  F  (9),  repre- 
senting, as  we  shall  see,  respectively,  the  contribution  due  to  reflection  from 
the  surface  of  the  sphere,  and  the  contribution  arising  from  (p-l)  internal  re- 
flections, or  p  refractions  in  the  sphere.  From  the  present  point  of  view,  the 
expansion  equation  (lO)  has  two  advantages:  It  unfolds  the  stationary  points  of 
the  integrand  so  that  a  given  Integral  In  the  sum  contains  at  most  one  station- 


ary  pointy  and  it  permits  a  ready  identification  of  the  multitude  of  terms 
in  acco27dance  with  ray  theory.  However^  the  integrals  are  not  yet  in  a  proper 
form  for  evaluation  because  we  must  distinguish  the  contributions  of  the  refrac- 
ted rays  and  the  diffracted  rays.   In  order  to  make  the  separation^  we  first 
extend  the  integrals  along  the  positive  real  v-axls  to  the  entire  real  v-axls. 
To  this  end,  we  observe  that 


H  -'   =  e^V^) 


(1) 

■V 


(2)   ^  e-^^'^H^^) 

■V  V 


-1/2 


-1/2 


^12,21^-^)  =  \2,21^^^ 


(16) 


From  these  illations,    it   follows   directly  that 


^  u 

/     dv   f         P     ,  /_(cos  9)   =      /      dv  f  -,    P     .  /o(cos   9)        for 

J  p,m     v-1/2^  J  p,p-m-l     v-1/2^ 


p  even. 


(17) 


The  terms  with  p  odd  have  a  slightly  different  reflection  property  about  the 
origin  of  the  v-axls.      To  reveal  It^    substitute  for  P        /   ( cos   9)   which^    from 
equation   (84)   of  the  appendix^    is  given  by 

^-1/2^'^°^   Q)    =   -  i^  '"''   Pv-l/2(-^°"   9)   +   (1  ;  e   ^2«^)    <^[]_}^i9)    •  (i8a) 


The  function  Cr  '(9)  is  defined  in  equation  (82).   It  is  a  linear  combination 
of  generalized  Legendre  functions  of  the  first  and  second  kind  and  is  designed 
to  have  a  purely  exponential  asymptotic  form.   In  view  of  the  sum  over  all 
values  of  m  that  occurs  in  F  (O),  the' second  term  on  the  right  yields  no  con- 
tribution when  the  substitution  for  P  _-i/Ji'^°^   ®)  ^^  ^^<^^  ^^  ^  (®^"  ^^'^^   "^^^ 
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remaining  term  in  (18a)  replacing  P^_^/2(co6  9),  it  is  clear  that 
00  o 


Combining  equations  (17)  and  (19)  with  (l^),  and  letting  the  summation  over 
m  commence  with  the  value  m  =  0,  ve  obtain 

p+1   «)    » 
^p(e)  =  -  ^  (-)  '  ^   f    ^^   ^p,m  '''""   ^^-1/2^--  «)'P  =  1^5,5,...  ,  (20) 


V^^  =  -  I  (-)'  ^  /  ^^  ^p.m  ^'""^  ^-1/2^-^  ^)>   P  =  2.^.6,...  .    (21) 


In  the  same  manner,  we  find  f or  P  (9)  that 


^o(^)  =  -  H  i  S  (-)"  /  ^-v  e^2«-  e^-  Vl/2(--  ^) 

In  order  to  separate  off  the  contribution  of  the  geometric  optics  fields 
we  make  use  of  an  idea  introduced  by  Franz  LJ_   Thus^  we  again  utilize  equa- 
tion (i8a)  to  substitute  for  P  i/p(<^os  G)  or  P    ,  (-cos  G)  into  F(G).   In 
addition^  the  following  alternative  to  (i8a)  will  also  be  used  which  follows 
from  (85)  of  the  appendix. 

P^_^/2(=°^  ®)  =  ^^"^''''  Pv-l/2(-=°"  e)  +  (1  +  e-^^")    <^[l}^{9)    .         (IBB) 


10  - 


When  either  of  the  substitutions  given  by  equations  (i8a),  (18b)  is  made  in 
F  (O)^  it  is  found  that  only  one  term  of  the  factor  with  Q  /  (9)  survives 
the  summation  over  m,  contributed  by  the  term  m  =  0: 

V.)=i(-)'^(/avf^^^e^(P-)-(.)ei^(^)(.) 


;  2]   r  ^v  f    e^^^   e^^'^^  P    ,Jcos 


P=  1,3,5....  ,  (23) 


^+1 


p-  2,h,6,...    ,  (2lt) 

where  the  choices  of  sign  are  not  independent,  all  upper  signs  belonging  to 
the  superscript  (l)  on  Q   /p(®)'  ^^^  ^^^   lower  signs  belonging  to  the  super- 
script (2).   The  upper  sign  corresponds  to  the  use  of  equation  (i8A)  and  will 
be  referred  to  as  choice  A.   The  lower  sign  corresponds  to  the  use  of  equation 
(i8b)  and  will  be  referred  to  as  choice  B. 

In  treating  the  reflected  field  In  this  manner,  only  choice  A  is  utilized. 
There  results 


F  (0)  =  F°(9)  +  F^(e)  ,  (25) 


h(')(x) 


'>)--i  /--niir-  ^^>)^  ^^'^ 


H.  '(x) 


1^(«)  =  i  2  (-)"  /  clv.v  e^^^-  e*"  P,.,/,(-=o.  9) 


The  first  Integral,  involving  0  i/o^^^'    ^^   called  the  geometric  optics 
field  and  is  denoted  by  the  superscript  G.   The  remaining  terms  are  called 
the  diffracted  field  and  are  denoted  by  the  superscript  D.   The  justification 
for  this  terminology  depends  on  the  evaluation  of  the  Integrals  and  ttieir 
geometric-optical  identification.   This  will  be  accomplished  in  the  next  sec- 
tion. We  do  not  yet  label  the  geometric  and  diffracted  field  contributions 
for  F  (0),  because  it  turns  out  that  firstly,  both  choices  A  and  B  must  be 
made,  depending  on  the  numerical  values  of  p  and  N.   Secondly,  it  is  necessary 
to  make  an  s-fold  substitution  for  Q  i/o^^^    ^^°^   equation  (i8a)  and  (i8b), 
again,  depending  on  the  values  of  p  and  N: 

s-1 

X  Z]  (-)Vi2rtev  ^   ^^  1,2,5,...  ,  (284) 

t=o 


^.%M  -   (-)==-'^=\i0.(«)*[f,.V^(=cse)-i,-^«\.,/,(-co.e)] 


X  2(-)V^='"^ 


l=o 


=  1.2,5, 


Finally,  then. 


Fp(e)   =  F^(0)   +  5^(9) 


:(«)  =  ^(-)^*7-^p,o'^'""""Vi)a<«> 


f"(9)  .  y-i 


/2 


(-COS9) 


(28b) 


(29) 


(50A) 


llt=0     ^  ' 

—  00 


choice  A 
P  =  1,5.5. 


»)  =  !(-)  = 


ig/-^<--'«%,.V./a(--' 


(51A) 


(30B) 


-S/-^<-"^^'-^„.p..,/a<--> 


choice  B 
P  =  1.3.5. 


(51B) 


15  - 


.>)  =  i(-)''°"/-^'--"Vi}.(«>^p,o 


2     Ave^ 


(p-2s)rtv 


^p,^^-l/2(=°^®^ 


P   ,  f     S-1         00 


^(e)  =  |(-)^       Z  /-^<-">"\,.p..v.(-: 


2  /-^(p--)-.^^^.^.^^^(_«: 


(52A) 


choice  A 


(55A) 


(52B) 


choice  B 
P  =   2,i+,6, 


(55B) 


We  see  from  equations   (l2),    (13),    (29),    and  (50)    -   (55)^    that  the  scattering 
amplitude  is   conveniently  decomposed  into  a  sum    of    fields  F  (9).      Each  of 
these  partial  amplitudes   is   self-contained  and  is   characterized  by  p,    the  nimi- 
ber  of  refractions  that  a  ray  has  undergone  before   contributing  to  the   scat- 
tered field  at  the   angle  e.      It  is   also  the  number  of  straight  line  segments 
of  a  ray  in  the  interior  of  the   sphere.      Each  peirtial  field  consists  of  tvo 


ll+ 


parts ^  a  geometric  optics  contribution  and  a  diffractive  contribution.   It 
suffices  to  know  how  to  evaluate  a  typical  partial  field  F  (0)  to  know  how  to 
evaluate  all  of  them^  although  the  reflected  field  F  (G)  is  a  little  differ- 
ent from  the  others.   We  are  now  in  a  position  to  evaluate  the  integrals  appear- 
ing in  equations  (50)  -  (35). 

IV.   The  Geometric  Optics  Field  for  the  Reflected  Ray 

To  evaluate  a  typical  Integral  appearing  in  F  (Q)^  we  replace  the  Hankel 
functions  appearing  in  the  integrand  by  their  Debye  asymptotic  form,  and  the 
0   /q(©)  functions  by  their  Laplace  asymptotic  form,  which  are  listed  in  the 
appendix,  equations  (88)  -  (91).  After  the  above  substitutions  are  made,  the 
integral  is  evaluated  asymptotically  by  the  method  of  steepest  descents.   The 
criterion  for  the  proper  separation  of  the  geometric  optics  field  from  F(0)  is 
just  the  existence  of  a  saddle  point.   In  this  manner,  the  choice  between  A 
and  B  and  the  value  of  s  can  be  detennined.   Only  the  major  contribution  of 
the  saddle  point  evaluation  will  be  retained.  ' 

Parenthetically,  we  remark  that  the  use  of  the  Debye  asymptotic  formulas 
always  appears  to  be  a  good  initial  step  in  treating  these  types  of  integrals 
involving  Hankel  functions.   Usually  it  is  true  that  the  major  contribution 
of  the  integral  does  come  from  a  region  whei^  the  Debye  approximation  is  valid. 

Using  equations  (88)  and  (90)  -  (93)  of  the  appendix,  we  obtain  from 
equation  (26) 

fQ(q)  ^  _iS  ^       /  sin  T  (cos  t)^  Re   °  dT,   0<9<)t,    (5^) 
in  which  we  have  observed  that 
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N  sin  T^-Bin  t   sin(T  -  t) 
N  sin  t  +sin  t   sin(T  +  t) 


Rll~H  = =  — ,  (55) 


0  =  2[]t  -  9/2)  cos  T  -  Bin  :] 
T  =  cos"  (v/x),   0  ^  V  <  X 
T  =  cos"  (v/Nx) 


(36) 


The  quantity  R  is  the  Fresnel  reflection  coefficient  appropriate  to  the 
tangent  plane  at  the  point  vhere  the  incident  ray  hits  the  sphere.   See  Figure  1, 
which  also  illustrates  the  geometrical  interpretation  of  v  as  an  "impact  para- 
meter" measuring  the  distance  of  closest  approach,  in  units  of  l/k,  of  the  ray 
to  the  scattering  center.   The  expression  "impact  parameter"  is  common  to 
quantum  mechanics,  in  which  the  incident  plane  vave  is  thought  of  classically  as 
a  stream  of  uniformly  spaced  particles  moving  orthogonally  to  the  surfaces  of 
constant  phase  of  the  plane  vave.   From  equation  (56),  we  see  that  the  restric- 
tions on  V  imply  that  we  consider  here  only  those  rays  which  impinge  on  the 
sphere  (fig.  l) . 

The  appearance  of  the  large  parameter  x  multiplying  the  phase  function  in 
the  integrand  in  equation  (5^+)  is  decisive  in  making  the  integrals  a  candidate 
for  evaluation  by  the  method  of  steepest  descents  (see  the  appendix).   The 
condition  for  a  saddle  point  is  the  vanishing  of  d0  /dT ,  or 


--°  =  -  2  sin  t(t  -  |)  =  0  ,  (57) 

dT  ■^ 

which  yields  the  saddle  point  value  (designated  with  subscript  o) 
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REFLECTED  RAY 


INCIDENT 
RAY      — ► 


TRANSMITTED  RAY 


FIGURE   I 


The  incident  ray  hits  the  surface  of  the  sphere  and  gives 
rise  to  a  reflected  ray  and  a  transmitted  ray.   The  radius  of 
the  sphere  equals  x,  in  units  of  the  wave  length  =  2jt/k.   The 
index  v  measures  the  distance  between  the  incident  ray  and  the 
center  of  the  sphere.  At  the  impact  point,  the  angles  t  and 
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,^  .  I  ,  ,^  .  ,  ecs  I  .  (58) 

The  solution  x  =  0  vhlch  arises  from  the  vanishing  of  the  factor  sin  t  must 
be  excluded  because  the  Debye  asymptotic  forms  lose  their  validity  there.   Fur- 
thermore, because  of  the  restrictions  of  the  Laplace  asymptotic  forms  for  Q  (O), 
the  forward  (9  =  0)  and  baclcward  (9  =  jt)  scattering  directions  must  also  be 
excluded.   These  directions  are  caustics  of  the  family  of  reflected  rays.   The 
saddle  point  is  of  second  order  (d  0  /dr  4-   O)^  and-  according  to  equation  (99)^  '"^ 
obtain 


^  \j^-   cos^(9/2)-  sin(9/2) 
\l^-   cos^(9/2)+  sin(9/2) 


(e)  --  -  ^  -——_—- : ,-2ixsin(9/2)^^  <  9  <  «  .   (39) 


This  is  the  geometric  optics  externally  reflected  field.      In  the   limit  in  vhlch 
N  ->  ioo  ,    this   reduces   to   the   reflected  field  of  an  Infinitely  conducting  sphere 
[j,l6],    as   it  should.      In  this   case,   the  phase   function  of  F     is  uneiffected  while 
the  reflection  coefficient  R       approaches   -1,    so  that  the  amplitude  of  F     is 
(-1/2). 


V.        Geometric  Optics  Field  for  the  Transmitted  Rays 

We  proceed  in  the  same  manner  for  the  p       transmitted  ray  as  given  by 


(50A,B)    and   (52A,B).      Then 


2,    sS+1 


F^O)^  JL-Llll-      r    sinT(cosT)2(l-R^)RP"^exp(ix0   (±9)±i  j)dT,  p=l,5,5,...    •    (^O) 
P  >/2nBin9    -^^  V     -^  / 

1  - 

i       s         2  1 

F°(9)-^^^y—     r    sinT(cosT)^(l-R^)R^"^exp   (lx0    [±(n-9)]±i  j)dT,  (i^l) 

^  v/2rtsin9    -^  V      P  y 


V^-    -      o 


p=2,l+,6,.-. 
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since  R   ~  R.   The  upper  and  lover  signs  in  the  exponent  go  together.   The 
upper  signs  correspond  to  choice  A,  vhile  the  lover  ones  correspond  to  choice 
B.  Although  the  limits  of  integration  in  equations  {ko)    and  (^l),  as  veil  as 
in  (5^),  extend  over  a  finite  interval^,  there  are  no  end-point  corrections  to 
the  saddle  point  evaluation  ve  are  maMng,  In  spite  of  the  restrictions  of 
equations  (56).   This  is  because  the  original  integrals  extended  over  the 
entire  v-axis.   The  phase  function  0  is  defined  by 

0  (±9)  =  2  JpNsinT^-sinx-  pi:^-T-(p-2s-l)|  +  I  costJ  ,  p  =  1,3,'^,...    ,  (^2) 

^    [>(:t-0)]=  2|pNslnT^-sinT-bT^-r-(p-2s-l)|  I  {n-Q)/AcoaA  ,    p  =  2,k,G, (^^5) 

The  condition  for  the  existence  of  a  saddle  point  for  these  phase  functions^ 
dd  /dx  =  0,  is  then 

9  =  +[T  n  +  2p(|  -  t^q)  +  2t^  -  2^    ,       V   =  1.5;5....  ,  (^^) 

g  =  ^  ±[7  „  +  2p(|  -  t^q)  +  2t^  -  2s7|  .   P  =  2>,6,...  ,  (i^5) 


inasmuch  as 


d0^(9) 


(6)  n  n   ©1 

2 =  2  sin  t  pT^-  T  -(p-2s-l)  2  "  2  • 


{^6) 


Assuming  the  saddle  point  is  of  second  order,  the  dominant  contribution  to 


fJ(9)  is 


Fj(e)  -  (-)I^"[Iine0^(T^j  %±n.^{^OB-:^)^^-\l-Aeyrp^0^{^^)±  l  +  ij 


p  =  1,2,3,...  ,       (^7) 
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vhere,  from  (^6), 


d  0  (t  )  /p  sin  T 

^-'(,  )=  -IEL^  =  2  sin  -:  (  i- ^  -  1   ,  (i.8) 

P  o      ^.2  o\   N  sm  T,^     '  ' 


0=0  (±9)        if  p  is  odd 


0  [±(n-9)J    if  p  is  even 


(^9) 


and  T  is  the  solution  of  the  transcendental  equation  (^^)  or  (^5)-   The  sign- 
in  front  of  the  factor  ^   in  the  exponent  in  (^7)  agrees  with  the  sign  in  ^-^i'^^)' 
It  has  been  assumed  above  that  the  amplitude  part  of  the  integrand  behaves 
like  a  constant  in  the  neighborhood  of  the  saddle  point.   Because  R  is  a  number 
which  is  always  less  than  unity,  the  larger  values  of  p  are  of  decreasing  im- 
portance because  of  the  quantity  R^"  which  appears.   The  latter  is  just  the 
attentuation  factor  due  to  the  repeated  internal  reflections. 

It  remains  to  specify  the  choice  of  sign  in  the  argument  of  0  as  well  as 
the  value  of  s.   The  criterion  for  this  is  Just  that  the  choice  should  yield  a 
saddle  point,  i.e.,  equations  (kk)   or  (^5)  should  yield  a  solution  0(t  )  such 
that  0  <  9  <  rt  and  0<t  -  ^-      It  is  easy  to  solve  the  saddle  point  equation 
explicitly  only  for  the  case  p  =  1.   Then  we  must  choose  8=0  and  selection  A, 
when  the  value  of  the  saddle  point  is 


.-1  N  sin(9/2) 

y/  -  I  cos(9/2)  +  1 


T  =  cos   .   \.  /     '  ^  (50) 


N  3ln(Q/2)        ^  ^^^^ 

^N^  -  I  cos(9/2)  +  1 
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From  (50)^  we  see  that  the  greatest  angle  of  refraction  occurs  for  t  =  0,  or 

o 

0  =  2  cos"-'-(i/n)  •  (52) 

This  angle  is  therefore  the  "botmdary  of  the  family  of  rays  associated  vith 
p  =  1  (see  Figure  2).   Because  of  its  association  with  t  =  0,  the  angle  at 
which  the  Dehye  formulas  become  invalid^  it  must  be  excluded  from  the  range  of 
applicability  of  the  partial  amplitude  F  (0).  From  equation  (^7),  the  contri- 
bution of  the  directly  transmitted  ray^  or  FHG),  is 


F^(9) 


N^-2Ncos(q/2)+1 
^1 


jjysec(e/2)  [Nco8(9/2)-i]  [n-cos(0/2)] 
2[n^-2Nco8(0/2)+i] 


X  expM.2x  v'N^-2Wcos(©/2)+l  j  ,  0  <  0  <  2cos""^  i  . 


(53) 


Except  for  the  factor  in  square  brackets  which  has  its  origin  in  the  transmis- 
sion loss  that  occurs  when  the  ray  enters  and  leaves  the  sphere,,  this  result 
agrees  with  the  differential  cross  section  for  the  scattering  by  an  attractive 
square  well  as  calculated  according  to  the  trajectories  of  classical  mechanics 
[17J  .  Naturally,  the  phenomenon  of  partial  reflection  is  alien  to  classical 
mechanics  so  that  the  difference  of  the  factor  (l-R  )  must  occur. 

The  saddle  point  evaluations  for  F,,  and  F^  are  essentially  the  same  as 

R      1 

those  made  by  Bucerlus  '-  -'  in  his  approximate  treatment  of  the  scattering  of 
light  by  a  water  droplet.   For  arbitrary  p,  it  is  not  possible  to  exhibit  t 
as  an  explicit  function  of  0.   Nevertheless,  an  Important  property  of  equa- 
tions (h-k-)    and  (if5)  can  be  stated:  there  is  always  a  saddle  point  at  t  =  jt/2 
when  s  =  0.   This  value  for  the  saddle  point  (v  =  O)  corresponds  to  the  ray 
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REGION 


FIGURE  2 


The  directly  transmitted  rays,  associated  with  the  value 

p  =  1,  are  confined  to  the  Illuminated  region  0  <  0  <  2t 

,-11  10 

where  -r^^  =  cos   -.   The  extremum  values  occur,  from  equa- 
tion (kh),    for  T 
be  chosen  equal  to  zero. 


and  T  =  jt/2.   The  parameter  s  must 


which  hits  the    center  of  the   sphere.      For  this  value,    "^-.r.  =  Jc/2,    and  p  odd 
requires  G  =  0,   while  p  even  requires  9  =  n.      This  means   that,    for  odd  p, 
there   is   always   some  neighborhood  of  the   forward  direction  which  is  illu- 
minated.     In  geometric  optics,    it  is  usual  to   call  this  the  lit  region.      In 
the  present   context,    lit   region  means  the  region  of  accessibility  of  the 
family  of  rays   associated  with  a  given  value   of  p.      This  region  is  different 
for  different  values   of  p.      Similarly,    for  p  even,    there  is   always   some 
neighborhood  of  the  backward  direction  which  is   a  lit  region.      In  short,    the 
lit  region  for  odd-p  rays  is  essentially  forward,    and  for  even-p  is  essen- 
tially backwards. 

The  saddle  point  equations,    {hk)    and   (^5),    for  p  ^  2,    can  be  studied  in 
a  qualitative  manner  in  order  to  decide  the   choice  of  equations  A  or  B  and  the 

value  of  8  for  t     •     We  shall  use  the  notation   (s,a)    or  (s_,B)   to  denote   a 
P 

particular  selection.     Let  us  examine   first  the  selection   (s,A)    for  even  p. 

In  the  neighborhood  of  t     equal  to   zero  or  jt/2,   we   find  from  equation   (k^)  with 

the  upper  sign  that 


0  =  2p  sin 


■1  1 


2rtS    +   €    + 


e  =  rt   -  2rtS  +  2(1  -  |)€  +   ...  when  t^ 


(5M 


Furthermore, 


=   0 


.1    17:1 


(55) 


These  relations   are  sufficient  to  describe  the  quaJLltatlve  behavior  of  9  as  a 
function  of  -r     for  various  values  of  N  and  p.      Th\is,    for  N  >  p,    9  is   a  monotonic 
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(2p  sin"     jj  -   2rts)    to  a  maximum  at  t    ,    and  then  decreases   to  the  value   (n   -   2jt8) 
at  T     =  rt/2.      This  behavior  is   illustrated  schematically  for   (0,A)   and  p  =  4  in 
Figure   3.      The   condition  that  9  be   restricted  to  0  <  0  <  n  nov  produces  the  fol- 
lowing selections.      For  N  >  p^    clearly  the   choice   (0,A)    is  to  be  made.      As   N 
decreases  below  p^    values   of  G  >  rt  would  appear  for  (0,A) .      Consequently  we  must 
choose  instead   (0,B).      This  has   the  effect  of  reflecting  the  e(T    )    curve   for 

(0,A)   in  the  line  9  =  -a,    and  is  suitable   for  esc  ~  ^  N  ^  v.      For  N  <  esc  ^, 

2p  2p 

9(t^)  for  (0,A)  becomes  larger  than  2n.   This  requires  the  selection  of  (l,A). 
As  can  be  seen  from  equations  (5^),  this  has  the  effect  of  reducing  9  by  2n. 
This  choice  applies  to  the  neighborhood  of  t  =0.   In  the  neighborhood  of 
T  =  n/2,   we  must  always  end  up  with  s  =■  0.   In  general,  for  p  large  and  N  close 
to  unity,  a  succession  of  selections  occur  as  t  varies  from  zero  to  it/2.  The 
sequence  commences  with  (s,A)  or  (s,b),  perhaps  increases  in  s-value,  but  ulti- 
mately decreases  and  concludes  with  (0,B).  A  typical  sequence  might  be 

j(8,A),  (3,B),  (8,A),  (s-1,b),  (8-1,A)...,  (1,A),  (0,B)|  .   In  Table  1,   we  sum- 
marize the  selections  as  a  function  of  N  for  p  =  2.   In  this  case,  we  must  take 
s  =  0  for  all  values  of  N.   The  lit  regions  are  also  indicated  in  the  table. 
For  p  =  ^,  the  selections  for  various  ranges  in  value  of  N  are  indicated  schema- 
tically in  Figures  31  -  3V.   These  figures  correspond  to  the  ciirves  labeled 
I  -  V  in  Figxire  3. 

The  consideration  of  equation  (hk)    for  odd  p  is  quite  similar.  We  have 
seen  that  for  p  =  1,  too,  the  selection  (0,A)  la  required  for  all  values  of  N. 
To  illustrate  the  situation  for  larger  values  of  p,  we  have  plotted  the  saddle 
point  equation  {hk)   in  Figure  4  in  a  schematic  manner,  for  p  =>  3^  s  =  0,  and 
selection  A,  for  a  few  values  of  N.   The  reduction  to  the  Interval  0  <  9  <  jr 
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TABLE   I 


N 

-^o 

IllTJiDinated 
or  lit  region 

Selection 

2  s  w  S  « 

°  <  ^o  ^  1 

h  sln"^  1  <  0  <  rt 

(0,A) 

■    ese  ?•  S  N  S   2 

°  <  ^o  ^  ^; 

i^-   sin'"''  i  <  9  <  rt 

(0,A) 

^;  ^  ^o  ^  1 

0(t^)    <  0   <  Jt 

(O.B) 

1  S  N  g   CSC  J 

°  <  ^o  ^  1 

©(Tg)    <  ©  <:t 

(0,B) 

Table  I .  The  selections  of  (0,A)  or  (0,B)  for  varying  values  of  N 
for  the  geometric  optics  field  for  p  =  2.   The  value  t'  is  that 
value  of  T  <  T_,  for  which  ©(t  )  equals  «,  as  given  by  equation  (k^) 

O      £5  O 

for  p  =  2  and  s  =  0.   The  angle  t^  =  sin"-'"[(N^-l)/5)]''"' ^  is  defined 
in  equation  (55) • 
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Tg=s.rn-'  [(N^-I)/I5]' 


FIGURE  3 


The  scattering  angle  9  as  a  function  of  t  as  determined 
by  equation  (^5)«   The  curves  are  schematic  and  not  literal, 
for  p  =  ^i  s  =  0,  and  selection  A  (upper  choice  of  sign).  A 
succession  of  values  of  N  are  shown,  varying  from  N  =  »  to 
N  =  1.   The  dotted  curve  labeled  r  is  the  locus  of  points 


for  which  dG/dx  =  0.   The  t  =  0  intercept  Is  9=8sln 


■1  1 
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e(To) 


27r  - 


-37r 


a^sin-i  [(n'-I)/8]^ 


FIGURE   4 


Schematic  Illustration  of  equation  (kk)    for  p  =  3^ 
8=0,  selection  A  (upper  choide  of  sign),  and  a  suc- 
cession of  values  of  N  from  N  =  oo  to  N  =  1.   The  dot- 
ted curve  Tg  is  the  locus  of  points  for  vhlch  de/dt  =  0. 

6  sin"-'-(l/N). 


t  § 
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via  the   correct  sequence  of  selections   Is   then  Indicated  In  Figures  4l   -  ^IV. 
The  values  of  N  for  which  each  selection  Is   appropriate  are  also  given.     As 
soon  as   a  value   of  s   Is  utilized  vhlch  Is  g  1,    the  lit  region  extends   over 
all  possible   scattering  angles.      The   geometric   significance  of  the  parameter 
s   Is   Jvist  this:   by  means   of  p  refractions   in  the  interior  of  the  sphere,    the 
p-ray  is   able  to  encircle  the  origin  s   times.      The  total  angle   of  refraction 
is  then,    so  to  speak,    9  +   2jt3.      This  is   the  reason  for  the  utilization  of 
equations   (28a)    and  (28b).      For  p  =   5  and  p  =  ^,    the  largest  value  of  s  that 
is  needed  is   s  =  1.      It  is   almost  superfluoixs   to  add  that  the   appropriate   form 
of  equation   (hf)  must   alvays  be  utilized  to  obtain  F   (9). 

VI.      Bovs 

It  sometimes  may  happen  that  the  saddle  point  is  of  third  order,  i.e., 
both  ^'  and  0"  vanish  at  some  point.   This  causes  a  greater  intensity  to  be 
produced  at  certain  special  angles,  and  Is  the  explanation  of  the  familiar 
ralnbov  that  occurs  in  the  scattering  of  light  by  water  droplets.   In  analogy 
with  this  phenomenon,  we  will  call  the  geometric  fields  associated  with  third 
order  saddle  points  bows.   If  we  Impose  the  condition  that  0"  given  by  equa- 
tion (^8)  vanish,  we  see  that  the  saddle  point  does  not  now  occur  for  every 
value  of  T  in  the  interval  (0,jt/2),  but  rather  only  at  the  special  value  t^, 

O  D 

wbere 


.  -1  M  -   1 


> 


P  >  1  (56) 


V      p2-i 
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The  bov  does  not  occur  for  p  =  1.   By  conrparison  vlth  equation  (55),  we 
see  that  the  bow  angle  9^  at  which  the  bow  occ\irs  is  always  a  scattering  angle 
of  local  minimum  or  local  maximum  deviation.   Therefore,  it  is  also  a  caustic 
for  the  set  of  p-rays.   Furthermore,  for  a  given  value  of  W,  the  bow  occurs 
only  for  those  values  of  p  for  which 


p  S  N  . 


(57) 


The  most  intense  bow  is  naturally  the  one  with  the  smallest  value  of  p  which 

satisfies  the  above  relation,  because  of  the  attenuation  factor  R   .   The  bow 

angle  is  obtained  from  (^4)  or  (^5)  with  t  =  t^.  At  the  bow  angle,  it  fol- 

O     B 

lows  from  equation  (lOO)  of  the  appendix  that  the  following  expression  appro- 
priate to  third  order  saddle  points  replaces  (^7)« 

1  1 


11  i 

F^CG^)  ~(-)^"x^(6Wn9g)    2r(i)  <JsinT(cosT)2RP-l(l-R2)[;0«(,)] 


(JE^p^^b)^-^) 


P  =  2,5.^,.. 


(58) 


where 


;„,2p,l„-l^iE!rJ)iJ\2sln-^  f€p^ 


1 


V   -1, 


.,,2psin-H  ff^]il  l-2^^-"'f  ^  )  -2- 


I^-^2 


P  =  5.5.7. 


,    p  =  2,1^,6,.., 


(59) 


■  -2k 


(t  ,9  )  =  2{/-lf{N^-lf   ,  p  =  2,^,k,...    ,(60) 


p'  B'  B' 


2  11 

2p(N  -l)8inT^cosT^    2  ^        2  2  2  2 


V^VV  =  5   3         =  ^(P  -1)  (P  -N")   ,  P  =  2,3,k,...    ^(61) 

N  sin  T^        p 

^^^b)  =FTT  '  P=  2,5,i^....  .(62) 


The  signs  in  front  of  ^  in  (58)  agree  vlth  the  signs  in  (59).   In  contrast  vith 
F  (9)  as  given  "by  (^7),  vhlch  is  of  order  \mity,  the  bov  term  as  given  by  (58) 
is  of  order  x  '  .   This  is  the  reason  for  the  greater  intensity  at  the  bov 
angle.   Naturally,  the  field  changes  gradually  from  order  x  '      at  the  bow  angle 
to  order  iinity  somevhat  away  from  the  bow  angle  on  the  lit  side.   To  understand 
this  transition  better,  we  must  examine  the  neighborhood  of  the  bow  caiistic 
more  carefully.   To  do  this,  we  expand  the  phase  fimction  0  in  a  Taylor  series 


1  ^\  5 

OT 

jre  — =2  is  given  by  (6I)  and 


cos  T^  P  =  5^5;,7^... 

cos  T^  p  =  2,l+,6, ...  , 

with  upper  sign  for  choice  A  and  lower  sign  for  choice  B.   The  first  improvement 
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in  F  (9)  in  the  neighborhood  of  9^^  comes  from  retaining  an  additional  term 
P  a 

in  the  exponent  of  the  integrand  of  the  integral  which  led  to  (58).   In  other 
words  we  retain  all  the  terms  in  the  Taylor  expansion  of  0  listed  in  (63) 
hut  no  more.   Then 

1  _  1  1 

F^(9)  ~  (-)P-'^x2(2nsln9)"  ^sinT^(  cost^)2rP-\t^)  [T-r2(x^7] 


X  exp       <ix 


S0, 


^w^^^w^«-v*f 


/ 


X     /      exp       Six 


5zf  '« ( T  ) 


sosr^  V  'b 


(65) 


The  integral  above  is  Alry's   integral,    equal  to 


x0;( 


6x 


1      ? 
2      \5    h\ 


V'^^b) 


^9^^%'' I 


(66) 


where  A  is  Alry's   function,    defined  by 


A(z) 


I  /  ^^[^i(y^  -  ^y}} 


dy   . 


(67) 


It  derives   its  name   from  Airy's  investigation  of  the  raLnljow  problem  over 
100  years   ago    L^^J .      Substituting  equations   (56) , (59-62),    (64)    and  (66)   into 
(65),   we  obtain 
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11  1 


n  2 


A9)  ^  (-)^^^(2nsin9)   \-3h^-lfiv^-^^)^h/-l)    ^Vi^) 


lH-1' 


'M 


exp        < Ix 


2(P^-1) 


11  i  _  i  -] 

2(/_i)2^(p2_jj2^2^p2_^-  2(Q_e^),  n\ 


{5p^(p'-l)-^(/-N^)}^]    , 


exp       < ix 


h(e-0^)x^  J  5p  (p  -1)     (p 


-11  i         _i 


P  =  5,5,7,... 


B^*  ¥ 


2  1^ 

:(9-9g)x5|5p^(p2-i)-2(p2_/)J5  ^ 


P  =   2>,6,., 


(68) 


vlth  upper  signs  for  choice  A  and  lower  signs  for  choice  B.     Equation  (68)   reduce* 
to   (58)    at  e  =  0_,,    as  it  should,    in  view  of  the  fact  that 


1   ^   „.l. 


A(z)'^-    -J-  r(^) 


(69) 


In  the  neighborhood  of  z  =  0.  For  |z|  »  1,  on  the  other  hand,  the  Airy  function 
decays  exponentially  on  the  shadow  side,  z  negative,  and  for  z  positive  behaves 
like 


A(.)-(f)^(|)'^=OB[^(f)=-|]  . 


(70) 
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Therefore  F  (0)  loses  its  x    dependence  and  behaves  like  a  constant  on  the 

2/5 
illuminated  side  of  the  bov.   This  occurs  vhen,  essentially,  x  '  (9-9^)  be- 

B 

comes  large. 

It  is  pertinent  to  inquire  whether  any  "super-bovs"    exists,    that  is  to 
say,   whether  there  exists   a  saddle  point  of  fourth  order.      This  would  require, 
in  addition  to  the  bow  saddle  point   conditions,    the  vanishing  of  0  '" .      From 
equation   (6I),    it  is  seen  that  the  possibility  exists  only  for  t     =  it/2.      There- 
fore,   from  (56)    and   {h-k)   or   (^5),   we  must  have 

N  =  p  and  3=0,      9=0,      for  p  odd  | 

f  (71) 

9     =  rt   ,      for  p  even 

The  existence  of  this  stationary  point  was   observed  by  Bucerlus    ^    -^ .      These 
directions  are   caustic  directions,    and  they  occur  at   angles   for  which  the 
Laplace   asymptotic  forms   for  Q  _-,/2(®)    ^-^"^  ^°'^  appropriate.      Consequently,    the 
super-bow  field,    just  as   the  field  in  the  neighborhood  of  a  caustic,    requires 
a  more  thorough  investigation,   which  we  shall  not  give  here.      We  shall  confine 

ourselves  to  the  observation  that   a  Debye   saddle  point  evaluation  of  a  fourth 

lA 
order  saddle  point  would  yield  for  x-dependence  of  the  super-bow  the  order  x 

It  may  be    conjectured    that  this   is  the   correct  x-dependence  becaiise  experi- 
ence with  a  similar   (inappropriate)  evaluation  of  the   forward  scattering  from 
Impenetrable  spheres  indicates  that  the  x-dei)endence  is   correctly  given,    al- 
though the   coefficient  is  not. 
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VTI.  The  Diffracted  Field 

We  will  now  evaluate  the  contribution  of  those  Integrals  which  we  have 
labeled  the  diffracted  field.   From  the  standpoint  of  a  stationaiy  phase 
evaluation,  these  Integrals  have  no  stationary  points.   Therefore,  it  sho\ild 
not  be  surprising  that  it  tiims  out  that  these  terms  are  exponentially  small 
compared  to  the  geometric  optics  field. 

It  is  to  be  kept  in  mind  that  the  diffracted  field  as  given  by  F^(9)  is 
applicable  only  to  those  regions  where  the  corresponding  field  F  (©)  is  valid, 
that  is  to  say,  the  illuminated  regions.  Although  the  evaluation  of  the  dif- 
fracted field  for  the  shadow  regions  would  be  essentially  unchanged,  addition- 
al contributions  to  the  diffracted  field  would  arise  from  those  Integrals 
labeled  geometric  optics  field. 

The  method  of  evaluation  of  a  typical  Integral  occurring  in  the  diffracted 
field  is  to  deform  the  path  of  integration  along  the  entire  v-axls  into  an 
infinite  semi-circle  In  the  upper  half  of  the  complex  v-plane.  The  factors  of 
e    in  the  integrand  are  sufficient  to  inspire  the  vanishing  of  the  integral 
along  the  semi-circle.   Therefore,  the  Integral  is  equal  to  2nl  times  the  sum 
of  the  residues  of  the  Integrand  lying  In  the  upper  halp  of  the  v-plane.   These 
residues  have  been  Investigated  by  Beckmann  and  Franz  ^  J  ^^  their  work  on  the 
dielectric  cylinder.  The  residues  are  due  only  to  the  zeros  of  the  function  L 
defined  in  equation  (ll).   These  zeros,  lying  in  the  upper  half  of  the  v-plane, 
are  given  by  the  expression 

i   ii 

V,  =x+«,(|)5e  5+  ...   ^  (72) 


T""^  ^i'^' 


where  q  are  numbers  which  are  determined,  in  the  first  approximation,  by  the 
equation 
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e   ^  (-)    w   ^  =  1  N  sin  t, 
^x'   A(q  )  1 


(75) 


A  more  complete  discussion  may  be  found  In  the  paper  of  Beckmann  and  Franz  (t.c.)' 

Let  us  consider  first  F  (e)  e£,   given  by  (51A,B)  or  (55A,B).   Due  to  the 
amplitude  factor  T,pRpp"  Tp-,  which  appears  in  f   ,    there  is  a  (p+l)-fold  root 
of  the  denominator  at  v  =  v.-   Consequently,  the  residue  sum  for  F^(9)  is 


^c 


P+-1  .  ^ 

rtX^  j=l 


B-1 


^(x)H^^^(x)H^-^^(Nx)H^^^(Nx) 


h(^)(x) 


"^  JiT 


(x) 


]•'■ 


||;2]'e^<--'''\../.(--)=""^j(i^) 


p-1 


^p 


vV 


p-1 


I(2)(.^  r.(i) 


3^    Vdv^     H^^^(x)H^^^(x)H^^^(Nx)H^^^(Nx)      H^ 

VVVV  V 


x)  Lsp^J 


X     ei(P-2^)-  e^^-  ^.i/2(=-^)    -'"''^y  (|^)''"' 


P  =  1,5.5. 


(7^) 


30  - 


(Nx)H^^^(Nx] 


V^^(Nx)' 
_v__ 

-  V   ^ 


^l(p-2s)jtv 


^-1/2^^°^®) 


-p-1 


00 


vv 


P-I 


V^^w(2^/',,-|TT'^^)/'M^\Tr(2)/,r^>,    tt(1) 


h[2)(,)  r-.d) 


h;  ^(x)h^  ^(x)h^  ^(nx)h;;  '(Nx)  H^"^(x)  L^v  ^^''^ 


(Nx) 


J^> 


v-1/2^      ' 


\  ,aLs-p-i 


P  =  2>,6, 


(75) 


with  upper  signs  for  choice  A,  and  lover  signs  for  choice  B.   The  quantities 
V  and  L  are  defined  in  equations  (ll).   Of  course,  the  value  of  s  in  these 
equations  is  already  prescribed  by  the  associated  geometric  optics  field  F  (9). 

These  terms  have  been  labeled  "creeping  vaves"  L'J  or  "diffracted  rays"  L J 
because  they  are  just  vhat  -would  be  contributed  by  rays  which  encircle  the  sphere, 
decay  exponentially,  and  shed  rays  tangentlally  as  they  move  along.   This  can  be 
seen  by  examining  the  phase  factors  appearing  in  (7^)  and  (75)-  With  the  appro- 
priate asymptotic  formulas  inserted,  the  exponential  factors  in  the  summation 


over  m  are 


31 


exp  <ix(2pNsinT-j_-2slnT)   +  iv   (p-2s)jt+(t+2T-2pT  ±9+2jtm   V         ,   p=   1,5,5,...    , 


(76) 


exp  <ix(2pNsinT   -2sinT)+lv   (p-2s):t+Tt+2T-2pT^±(jt-9)+2jtm   >         ,p=  2,4,6,...    ,        (77) 

I  Hi 

vhere  the  choices  in  sign  are  Independent.  The  +jt  arises  from  the  A  or  B  selection. 
The  ±9  and  ±(rt-9)  are  due  to  the  tvo  possible  exponentials  in  the  asymptotic  forms 
of  P    ;p(cos9)  and  P    /  (-cos9),  respectively.   Similar  expressions  are  contained 
in  the  sxzmmation  over  I.      From  equation  (72),  it  can  be  seen  that  at  v  =  v.,  t  '^-' 0 
and  T  '^  cos"  — .   Consequently,  the  diffracted  field  contribution  arises  from  the 

grazing  incident  ray  as  was  first  pointed  out  by  Keller  I 1  in  the  application  of 

geometric  diffraction  theory  to  infinitely  conducting  obstacles.   The  first  term  in 
the  brackets  describes  the  phase  factor  caused  by  the  p  transmissions  through  the 
interior  of  the  sphere.  From  the  second  term,  we  see  that  these  terms  are  exponen- 
tially small,  becaiiBe  v.  has  a  real  imaginary  part.   The  exponential  decay  factors 


exp  i-q  (|)5coB  Ijapsin'-'-  i+n-23:t±9+  2rtmJ  I       , 
exp  |-qj(g)5cos  |f2psin"^  i  ^  ^   _  2sn   ±   («-9)   +  2«^  I 


P  =  1.5,5,.-.  ,  (78) 


p  =  2,1^,6,...  ,       (79) 


The  quantity  in  the  square  brackets  is  precisely  that  angle  of  encirclement  about 
the  sphere  necessary  to  permit  a  tangential  ray  to  be  shed  in  the  9-dlrection, 
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where  0  Is  an  angle  in  the  lit  region.   It  is  to  be  kept  in  mind  that  the 
parameter  s  vas  chosen  so  as  to  insure  this.   In  contract,  the  quantity  2-im. 
describes  the  additional  decay  due  to  m-fold  encirclement  of  the  sphere.   The 
±  signs  in  front  of  the  quantities  9  and  (n-0)  that  occur  in  (78)  and  (79), 
which  have  their  origin  in  the  two  exponentials  appearing  in  the  Laplace  asymp- 
totic forms  of  P    /  (±cose),  describe  the  two  possible  ways  for  diffracted 
rays  to  be  shed  in  the  0  direction.   This  is  Illustrated  in  figures  5  and  6 
for  m  =  0,  p  =  1,2,  and  s  =  0. 

The  diffracted  field  F^(0)  which  is  associated  with  the  reflected  ray 
contains  a  sum  which  is  independent  of  R  ,  according  to  (2?).   This  sum  may 
be  shown  to  vanish.   This  follows  because,  on  the  one  hand,  the  imaglnajry 
part  of  the  integral  in  the  summand  is  an  even  function  of  v  while  the  real 
part  is  odd.   Therefore,  the  integral  over  the  entire  v  axis  is  pure  imaginary 
and  equal  to  twice  the  value  of  the  integral  from  zero  to  infinity.   On  the 
other  hand,  the  latter  Integral  has  a  null  imaginary  part.   This  may  be  seen 
by  deforming  the  path  of  integration  to  the  positive  imaginary  v-axis.   In 
view  of  the  fact  that  P     ;  (-cosO)  is  a  real  function,  it  follows  that  the 
integral  is  a  real  function. 

Thei^fore,  the  diffracted  field  F^(0)  consists  only  of  the  second  sum  in 
(27).   This  sum  may  be  treated  in  the  same  manner  ae  F{Q).     We  find  for  it 
that 


^<«>=-s<->"g  (f  £s^^"w---^^""^ 


(80) 
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FIGURE  5 


The  diffracted  rays  which  are  shed  in  the  direction  9  of  the  lit  re- 
gion for  p  =  1  and  m  =  0  are  given  by  Fr(9),  equation  (7^)  with  the 
selection  (0,A).   There  are  two  different  possible  angles  of  encircle- 
ment of  the  sphere  which  yield  a  diffracted  ray  in  the  direction  0. 
They  correspond  to  the  two  possible  exponentials  appearing  in  the  asymp- 
totic form  of  P    /_(cos  9).   The  angles  of  encirclement  are  indicated 
by  equation  (78)  with  choice  A.   Only  the  upper  grazlng-incldent  ray  is 
shown. 
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DIFFRACTED 
RAY 
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LIT  REGION.  6=  4cos-«i. 


FIGURE  6 


The  diffracted  rays  shed  in  the  direction  9  of  the  lit 
region  associated  with  p  =  2  and  m  =  0  are  expressed  by  F?(g)^ 
equation  (75)'   The  figure  illustrates  the  selection  (0,A). 


to  be  made.   The  decay  factors  corresponding  to  the  figure  are 
given  by  equation  (79). 
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This  expression  applies  to  the  lit  region  for  the  reflected  field  vhich,  as 
we  have  seen,  extends  everywhere  except  the  two  caustic  directions  9=0  and 
9  =  rt.   It  contains  a  decay  factor  given  by  equation  (79)  with  p  =  s  =  0  and 
selection  A-   If  we  pass  to  the  limit  of  the  infinitely  conducting  sphere, 
equation  (8o)  reduces  to 


nt=o       J=l      dH^  '(x)  ' 


(81) 
The  zeros  v.  are  the  zeros  of  H^   (x)  in  this  case.   This  agrees  with  the  re- 


*- -*    and  T.ew  nnd  Kella-r  L.-'J 


suit  of  Franz  *-  -'  and  Levy  and  Keller 

VIII.  Axial  Caustics  and  Boundaries  of  Illuminated  Regions. 

The  results  we  have  obtained  are  not  valid  in  the  neighborhood  of 
the  axial  caustics,  that  is  to  say,  the  forward  and  backward  scattering 
directions,  and  in  the  neighborhood  of  the  boundaries  of  the  illuminated 
regions  for  each  of  the  partial  amplitudes  F  (9).   In  the  latter  case,  this 
is  because  the  boundary  is  determined  by  the  saddle  point  value  t  =0,  and 
the  Debye  asymptotic  formulas  are  not  valid  at  thet  point.   This  state  of 
affairs  can  in  fact  be  easily  corrected  with  the  use  of  the  uniform  asymp- 
totic formulas  for  Hankel  functions  in  terms  of  Airy  functions.   While  the 
breakdown  of  the  Debye  formulas  may  also  contribute  to  the  lack  of  validity 
of  our  resiats  In  the  neighborhood  of  the  axial  directions,  the  more  serious 
cause  is  the  infinite  behavior  of  the  generalized  Legendre  functions  q}  '^      (9) 
in  the  forward  and  backward  directions.   In  this  case,  it  appears  necessary  to 


-  "^  - 


rewrite  the  expressions  for  F  (9)  so  that  they  are  valid  in  the  neighborhood 
of  these  directions.   This  will  be  investigated  in  a  subsequent  paper. 
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Appendix 


Useful  Formulas 

In  connection  with  Legendre  functions,  it  Is  convenient  to  define  q}  ^(9) 
and  Q^  (9)  in  the  following  way  so  that  these  functions  have  purely  exponen- 
tial asymptotic  forms 

Q^^)(9)  =  I  [p^(cos  9)  +  i  I  Q^(cos  97|  ,  (82) 

q^2)(9)  =  I  fp^(co6  9)  -  i  f  Q^(cos  97]  ,  (85) 

Here,  Q  (cos  9)  Is  the  generalized  Legendre  function  of  the  second  kind.   Then 
the  relation  between  Legendre  functions  may  be  written  as 

P  (-  cos  9)  =  e"^'^^  P  (cos  9)  +  21  sin  jtv  Q^-^Hq)   ,  {^) 

P^(-  cos  9)  =  e^'^''  P^(co8  9)  -  21  sln^v  (^^^(9)  .  (85) 

[211 
The  Laplace  asymptotic  forms  for  Legendre  functions  ^    -■  are 

_1 
P^(cos  9)  '^  (i  jtv  sin  9)  ^  cos  (v  +  |)9  -  J  ^  0  <  9  <  j(  ,   (8^) 

_1 
I  Q^(cos  9)~  -(|  «v  sin  9)  ^  sin|7v  +  i)9  -  f]  ,  0  <  9  <  s  •   (8?) 

With  the  use  of  theae  expressions,  the  sBymptotic  forms  of  qJ,^^'^^^(9)  become 
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q[^\9)  ~  (2rtv   sin  0)    ^  expTiCv  +  i)9  +  1  |^  ,     0  <  0 


(88) 


_  1 

Q^^)(0)  ~   (2rtV   sin  e)    ^  exp[+i(v  +|)0-ij^,     O<0<Jt     .  (89) 

r22i 

The  Debye  asymptotic   forms   for  the  Hankel  functions    >-   -"    and  their  deri- 
vatives,   for  large  eirgument  x  and  index  v  «  x,    are 

_  1 
h(^)(x)  '^^(|  rtx  sin  t)   ^  exp^x(Bin  t  -  t   cos  t)    -  i  J L  (90) 

_  1 
H^^^(x)  ~  (|  rtx  sin  t)   ^  expTixCsin  -r  -  t   cos  t)   +  i  jJ  ,  (91) 

<3B     (x)  1  "9       r  71 

—11^^^ ^i  sin  ^(i  nx  sin  t)   '^  exp  ix(sin  t  -  t   cos  t)   -  i  J    ,  (92) 

dH(2)(x)  -i 

V      ^ 


dx 


^-1  sim(^  jtx  sinT)        exp  -lx(sin  t  -  -r   cos  t)   +  i  ^    ,  (93) 


vhere  t  =   cos'     —     •  (9^) 


According  to  the  method  of  steepest  descents  *-  <J  ^  if  the  Integral  -0.  to 
be  evaluated  asymptotically  for  \  »M  is  of  the  form 

^  =  jT  A(|)e^^^(^)d5  (95) 
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and  contains  an  Interior  point  along  its  path,  of  integration  such  that 


M=o.   ^-0,...    0W_-^^O     at 


d| 


I  =  I, 


(96) 


<i| 


df 


then  it  is  said  to  have  a  saddle  point  of  order  m  at  |  =  |  .   If  0  and  A  are 
expanded  in  Taylor  series  about  the  point  |  =  |  ,  e.g.^ 


=  0(i,)*^/"'(g(i-g"+... 


(97) 


and  A  is  of  the  fonn 


A  =  A^+A^(g  -  1^)  + 


(98) 


then  the  dominant  term  (qua  function  of  X)  in  the  asyniptotic  expansion  of  «ft  is 
given  by 


i^0(^o) 


1 


"   ^  Te' 

-i    1X!ZS(|  )  r   6   "15   1 


if  m  =  2 


(99) 


)       if  m  =  5   .      (100) 
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